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Abstract. 

We present a scheme by which projective homodyne measurement of a microwave 
resonator can be used to generate entanglement between two superconducting charge 
qubits coupled to this resonator. The non-interacting qubits are initialised in a product 
of their ground states, the resonator is initialised in a coherent field state, and the 
state of the system is allowed to evolve under a rotating wave Hamiltonian. Making 
a homodyne measurement on the resonator at a given time projects the qubits into 
an state of the form {\gg) + e~*'^|ee))/-\/2- This protocol can produce states with a 
fidelity as high as required, with a probability approaching 0.5. Although the system 
described is one that can be used to display revival in the qubit oscillations, we show 
that the entanglement procedure works at much shorter timescales. 
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1. Introduction 

Central to the construction of any useful Quantum Information device will be an array 
of quantum systems whose collective quantum state can be prepared, measured or 
otherwise manipulated [1]. For the case of two- level systems — qubits — and in the 
solid state arena, small superconducting grains coupled to bulk superconductors via 
Josephson junctions, often referred to as Cooper-pair boxes, have been shown to be 
promising candidates for playing the role of the qubits[2]. Quantum state control 
has been demonstrated in single qubit devices through coherent or Rabi oscillations, 
observed with readout schemes consisting of either a single electron transistor (SET) [2] 
or quantronium[3] circuit coupled to the qubit being measured. Corresponding single 
qubit quantum state control results have also been obtained in superconducting flux 
qubit devices[l]. 

More recently, the coupling of such superconducting qubit devices has also been 
achieved[5l El [7]. Such coupling generates entangled two-qubit states, which is both of 
vital importance to possible quantum computing applications [T| , and also of interest 
in testing the fundamental limits of quantum mechanics in macroscopic objects[8], [9]. 
The two charge qubits of Pashkin et a/. [5] were capacitively coupled, and each qubit 
independently measured with SET readouts. The resulting two-body density matrix, 
was consistent with the existence of non-classical correlations between the two qubits, 
although decoherence limited the fidelity of the CNOT gate which was attempted in 
the two qubit device [ID]. Entanglement between superconducting flux qubit devices 
was reported by Berkley et a/. [6]. A more recent experiment [7] showed antiphase 
oscillation of a two-qubit system, and detailed quantum state tomography established 
entanglement with 87% fidelity between two macroscopic superconducting charge-phase 
qubit devices pT]. Significant progress has also been made coupling superconducting 
charge [12], [13], [H] and flux [15] qubits to a quantum resonator mode, demonstrating 
effects such as the AC-Stark effect and resolving individual photon number states in 
the resonator. As well as being an essential part of any quantum computing protocol, 
the generation of entanglement is necessary for the demonstration of a Bells' inequality 
violation [16], which would prove unequivocally that the system cannot be described 
classically [T7]. 

In all these experiments the coupling of the superconducting qubit to its 
measurement device has been shown to be one of the central aspects controlling its 
decoherence. Very long relaxation times were observed in the experiments of Wallraff 
et fl/.[T2l [T3l [T4j in which the Cooper-pair box was coupled to a microwave stripline 
resonator, but not to any other measurement device. Measurement of a resonator 
dispersively coupled to a superconducting qubit has thus been shown to provide an 
excellent 'readout' system for the superconducting qubit [T3l ITS] . Quantum resonator 
modes therefore have real potential for both coupling and readout of superconducting 
qubits. 

In this paper we present a measurement protocol by which entanglement may be 
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generated between two such superconducting qubits, which are coupled only though 
a common stripline resonator. We show that projective homodyne measurement on 
the resonator field may perform a projective measurement on the coupled qubit/field 
system, resulting in the generation of entangled states of the two qubits. It is interesting 
that there is assumed to be no direct coupling between the qubits themselves, and 
it is the act of projective measurement on the field which creates the entangled 
state. Our scheme for generation of entanglement is therefore analogous to the Knill, 
Lafiame, Milburn or 'KLM' protocol[l9] to generate entanglement of photons in linear 
quantum optics through projective measurement. In the same way, our scheme uses 
measurement to generate entanglement between qubits which have no other direct qubit- 
qubit interaction. The analysis presented here will be relevant for other physical systems 
consisting of two-level systems coupled via an harmonic mode such as atom-chips or 
nanomechanical systems. 

Our work builds upon the extensive theoretical study of Meunier et al. [20] of collapse 
and revival phenomena of two level quantum systems coupled to a single quantized 
radiation mode. As is stressed by Meunier et a/. [20] the rich variety of quantum 
phenomena displayed by the 'one qubit one mode' system is well known from Quantum 
Optics where the role of the two level system is played by a two level (Rydberg) atom [21]. 
A particularly striking example of these is the collapse and subsequent revival of the 
initial Rabi oscillations in the atom (qubit) section of the Hilbert space when the initial 
radiation field is in a coherent state |a). Remarkably, as was first noticed by Gea- 
Banacloche [22], at times t, between the collapse time tc and revival time tr-, when 
there are no Rabi oscillations, it is the radiation field part of the wave-function which 
manifests complex quantum behaviour. Following these insights we have studied the 
time evolution of the two qubit system coupled to a radiation mode. We show that 
simple projective homodyne measurements on the radiation field, at times between 
collapse and revival of multi qubit Rabi oscillations, may be used to obtain interesting 
entangled states of the two qubits. 

The paper is structured as follows. In section 2 we introduce the two qubit version 
of the Jaynes-Cummings model [231 l2l] and, in the interest of clarity, recall the current 
state of understanding of the complex evolution with time of a combined qubit and 
field state. We also discuss the expected results of a homodyne measurement on the 
field variable of our system. The main results of this paper are developed in section 3 
where we describe and analyse a simple experimental protocol which involves homodyne 
measurements on the field variable and results in heralded maximally entangled states 
of the qubits. Our general conclusions are presented in section 4. 

2. Jaynes-Cummings model for two qubits coupled to a resonator 

The salient features of one qubit coupled to a single mode of the radiation field can be 
described approximately by the Jaynes-Cummings model [23], much studied in quantum 
optics. This model has generated a great deal of interest in the past as it exhibits 
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interesting behavioural] and it has been used to describe quantum correlation. 

In this paper we consider two charge qubits each coupled capacitively to a striphne 
resonator using the two qubit Jaynes-Cummings modelpl]- If the qubits are placed at 
an antinode of the fundamental harmonic mode of the resonator, then we can describe 
the system as a pair of two-level systems coupled to a simple harmonic oscillator. The 
charging energy of the qubits and their coupling to the resonator can be controlled by 
the application of magnetic and electric fields[2]. If these are tuned so that the qubits 
are close to resonance with the field we can describe the system using a rotating wave 
Hamiltonian, 



(1) 



where 



a^=\e){e\-\g){g\, a+ = \e){gl a- = \g){el (2) 

and d are the creation and annihilation operators of photons with frequency u, Ei^2 
are the charging energies of the qubits, and Ai,2 the resonator-qubit coupling terms. We 
consider the system to be exactly on resonance, for simplicity, and also assume that the 
Cooper pair boxes have equal charging energy, (-Ei,2 = E) and coupling to the resonator, 
(Ai,2 = A). In the rotating wave approximation, we notice that states with m excited 
qubits and n photons only couple to states with m' + n' = m + n. This means that we 
can describe the system as an infinite set of non-interacting four level sub-systems all 
labeled by -|- m. We can find the solution of each of these sub-systems analytically, 
and then sum over these to obtain the state of the whole system. 
For an initial state, 



(3) 



n=0 



we obtain, 



oo 



n=0 



^yn{n - 1) 
2n- 1 



cos At V2(2n - 1) - 1 |ee, n - 2) 



- t-^=0==sm (^Xty/2{2n - 1)) {\eg,n-l) + \ge,n-l)) 
+ (n cos (AtV2(2n - 1)) + {n - 1)) \gg, n) 



(4) 



where the state \eg,n) is the tensor product of a Fock state \n) of the photons and the 
state of the qubits which is either excited (e) or ground {g) for each qubit. For brevity 
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we only discuss the solution for the case where the qubits both start in their ground 
state, but we note that the solution with a general initial condition is generically similar. 
The initial state of the field is determined by the coefficients C„, which for a coherent 
state I a) are given by, 

= e-l"l^/^^, (5) 

with the phase, 6, and average occupation, n, of the coherent state determined by 
a = y/Ee~'^^. 

In spite of its relative simplicity equation (jl]) describes a wide range of interesting 
phenomena. A much studied example of these is the collapse and revival of Rabi 
oscillations. In this paper we shall be concerned with phenomena at times where there 
are no Rabi oscillations. To illustrate the generic behaviour of {ipit)) in equation (jll) we 
shall now comment on the qubit and the photon sector separately. If we consider the 
state of the qubits after tracing out the field state, the probability of the qubits being 
in the state \gg) is, 

oo 

P,,{t) = {gglp'^mg) = J^liaaMm)]" (6) 

n=0 

where p'^it) is the reduced density matrix at a certain time for the qubits when the field 
has been traced out. In agreement with Iqbal et a/. [25], we find the results depicted 
in figure [TJ Evidently we see that initial oscillations of the qubit states decay rapidly, 
after which the qubit is in a mixed state. After a period of time, the oscillations revive, 
indicating the information about the qubits' initial state has returned to the qubits from 
the field. We see a larger revival at the revival time, = 27ry/E/X, which is proceeded 
by a smaller revival at tr/2. The smaller revival is not seen for the one qubit case and 
this is due to the extra frequency introduced by the addition of the second qubit. 

So far we have studied the time evolution of the reduced density matrix of the 
qubits. To investigate the radiation field it is useful to calculate the Q function [26j . 

Qia,t) = {a\p^{t)\a) (7) 

where p^{t) is the reduced density matrix at a certain time for the radiation field when 
the qubits have been traced out. The values of Q{a,t) in the complex alpha plane at 
fixed times are shown in figure [2l 

The Q function shows a set of three 'blobs', each representing a mesoscopic 
wavepacket of the cavity field. As a function of time the three wavepackets move around 
the complex plane and follow the circular path of radius v^- Although the wavepackets 
all begin in the same place, they evolve with different frequencies, so the states begin to 
separate. After a period of time, depending on the differences between the frequencies, 
the different sub-distributions are separated by more than their diameter and can be 
easily distinguished, (figure [2]d). 

To elucidate the connection between the collapse and revival in the qubit part of 
the Hilbert space and the peaks in the Q{a, t) distribution depicted in figure [2] we recall 



Entanglement of superconducting charge quhits by homodyne measurement 



6 




Figure 1. Revival of the initial state of the qubit system. The curve shows the 
probability that the qubits are in the state \gg) after tracing out the resonator. The 
initial qubit oscillations rapidly decay, indicating the qubits are in a mixed state. The 
revival of the oscillations indicates that the information about the qubit states has 
been transferred back to the qubits from the resonator. The main revival occurs at 
tr — 27rv^/A, with a partial revival at tr/2. We have chosen w = A = 1 and n — 30. 



briefly the results of Meunier et a/. [201 obtained in the large n limit. Generalising the 
one qubit result of Gea-Banacloche for the two qubit case they found that {ipit)) can be 
written as a superposition of Gea-Banacloche states described below. 



1 

k=-l 

where k=-l, 0, or 1, \Dk{t)) is a state of the qubits and is the state of the field. 

\D-m = l(e"^^(--;^ V^«>e) - e-'^--V^^'e-^%\eg) + \ge)) + \gg)) 

(9) 

1^ o(t)) = lilgg) + e-2^(''+-/2+-*)|ee)) (10) 
\D ,{t)) = i(e"^^(-+^)*e-^^^|ee) + e^'^^^V^^^'e-^Weg) + \ge)) + \gg)) 

(11) 

= e-*'=^^*|e"'^"+'^^)*«) (12) 

The state of the field in each Gea-Banacloche state \Dk{t)) ® |$fc(t)) is one of three 
coherent field states |$fc(t)), each of which has a phase that evolves with a frequency 
given hy to + k-^. Each field state corresponds to a particular qubit state \Di:(t)), as 
given by equations (pi)- (fTTl) . In particular, we note that qubit state \D o) corresponding 
to the field state |$o(^)) is a superposition of \gg) and |ee) only, and has no components 
of the states \ge) or \eg), a fact that we shall exploit to produce our entanglement 
protocol. 

The Gea-Banacloche states can be clearly seen in figure [2], where the three 
'blobs' represent the states \^k(t)), and we see that the frequency at which these 
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Figure 2. Phase-space plots of the Q-function for four different times, where the 
dimensionless x = a + a* represents the electric field and y = {a — a*)/i represents 
the magnetic field. Also plotted are white diamonds corresponding to the evolution 
of the field states in equation (fT^ . (a) t — tt/{2uj): The phase space close to 
t = 0, when the peaks are not well separated, (b) t = 37r/(2a;): After a period of 
time the sub-distributions, evolving with different frequencies, separate and can easily 
be distinguished, (c) t — tr/2: The time at which there is the first occurrence of 
spontaneous revival, (d) t = tr'. Location of the states at the second main revival 
peak. For all the diagrams A = w = 1 and n = 200. 

sub-distributions move around the complex plane is determined by the frequencies 
in equation f|T2|) . The well known collapse and revival phenomena, seen in figure 
[H can now be understood in terms of the evolution of the Gea-Banacloche states. 
The Rabi oscillations collapse when the field states are well separated (figure EJ^b)). 
At the large revival peak (figure El^d)) all the distributions overlap in the complex 
plane, which occurs when all the field states are in phase. The time this occurs is 
determined by the condition exp{iujt) = exp(i(co' -|- X/y/E)t) = exp(i(co' — X/^/E)t). 
This condition is fulfilled by tr = 2tt\%/\. At an earlier time we observe a smaller 
revival peak, when two of the states are in phase (figure EJ^c)). This occurs when 
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exp{i{u! + X/y/E)t) = exp{i{uj — \/^/n)t), namely at t = tr/l. 

After gaining a clearer understanding of the physics by taking the large-n limit, we 
now return to our exact calculations based on equation (jl]). Whilst the Q-function shows 
a fuller understanding of the state at a specific time, we shall study an experimentally 
more accessible distribution defined as follows. Evidently the eigenvalues of the operator 
X = (a + a^) are related to the position variable (electric field) of the harmonic oscillator 
which describes the cavity mode. In order to analyze the behaviour of this system as a 
function of time, we can plot the corresponding a;-distribution of the field state after we 
have projected out the qubit states, i.e. we plot 



where r = gg,ee,eg,ge (figures [3t^a) - [3](d)). A convenient numerical method for 
obtaining this distribution is to project out a particular qubit state, transform the 
density matrix into the x-basis and take the diagonal elements. Obviously, the finite 
number of basis states in the calculation will lead to a discrete set of position eigenvalues 
Xi for X, and hence a discrete set of probabilities Pr{xi,t), but this set of points can be 
made sufficiently dense to form a probability density, as given in equation f|T3l) . 

Although the Q{a,t) function studied by Meunier et al. [20] gives a fuller account 
of the cavity field than the above Pr{xi,t), nevertheless, as shown in figure [3], the latter 
also captures the salient features of a very interesting quantum state at hand. In the 
interest of clarity the pictures in figure [3] are at the same times as the pictures in figure 



As in figures [2]^a) -[2](d) we can see that although the distributions initially overlap, 
they evolve with different frequencies, so the peaks begin to separate (figure [3](b)). After 
a period of time the distributions overlap and cause the two revival peaks (figures [3]^c) 
andEKd)) 

3. Entanglement protocol 

The basis behind the proposed entanglement procedure is a measurement of the quantum 
field at some chosen time, which projects the state of the two qubits into an interesting 
and useful entangled state, heralded by the outcome of the field measurement. In 
quantum optics the idea of measuring one part of an entangled system to learn something 
about the other goes back over twenty years [271 EH] • Measurement of a coherent state 
entangled with a photon can be used to project the photon state without absorption 
of the photon. In our case there is a tripartite entangled system and projective 
measurement of the coherent state is used to project the remaining bipartite qubit 
system. The relevant measurement is balanced homodyne detection, where the coherent 
state is mixed with a strong local oscillator field on a 50:50 beam-splitter before photo- 
detection of the two outgoing fields. It is well known, e.g. [29], that in the strong 
local oscillator limit the photon counting measurements made on the two outgoing 
modes of the beam-splitter correspond to a projector (e.g. onto a chosen 
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Figure 3. The x-quadrature distribution of tlie field state after projecting out 
the qubits states r = \gg),\ee),\ge) + \eg)/\/2 for a series of different times. The 
vertical lines indicate the location of the field states as given in equation ([T2|l . (a) 
t = 7r/(2a;): The x-distribution close to t = 0, when the peaks are not well separated, 
(b) t = 3tt/(2uj): After a period of time the sub-distributions, evolving with different 
frequencies, separate and can easily be distinguished, (c) t — tr/2: The x-distribution 
at the first revival peak, (d) t — 1^: The x-distribution at the second main revival peak. 
The vertical lines indicate xit) for harmonic oscillators with frequencies w, w ± Xj^ffi. 
For all the diagrams A = a; = 1 and n — 200. 



quadrature of the initial coherent state field being measured, where the quadrature 
is set by the chosen phase of the local oscillator field. As will be seen for our system, 
at certain times a projection onto the x quadrature of the field entangled with the 
qubits can leave the qubits in an entangled state. For the case of superconducting 
qubits coupled to a microwave field mode [30l |31] this requires homodyne detection 
at microwave frequencies, rather than the familiar beam-splitter and photo-detection 
systems employed at optical frequencies. The mixing (of signal and local oscillator) 
and detection at microwave frequencies, in order to project onto a defined quadrature 
of a microwave coherent state field, may be achievable through use of a single electron 
transistor (SET), as has been discussed in detail in [32] . 
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We stress that in our entanglement protocol, the idea is to introduce, or turn on, 
the homodyne measurement applied to the cavity field at some chosen time. This would 
require some sort of fast (on the time-scale of the qubit /field evolution discussed in 
the previous section) opening or gating of the cavity, to enable the microwave field 
to be subject to measurement at a chosen time. Given the long coherence times (or, 
correspondingly, high quality factors) observed in superconducting qubit and microwave 
mode experiments |12l [13], HI] , we neglect decoherence in the field mode and qubits prior 
to the measurement, but explicitly allow for the decoherence of the actual measurement 
process, through its projection of the state of the qubits and field. 

So, taking homodyne measurement to act as a projective measurement of x, our 
entanglement procedure is as follows. At a given time t the x-quadrature measurement 
will give a value of x taken from the distribution P{x,t) = Pr{x,t). Repeated 



observations would measure the whole distribution P{x, t). Now if the peaks in Pr{x, t) 
are well separated for different qubit states r, then the result of a single x- measurement 
indicates which sub-distribution Pr{x,t) the system is in, and thus the value of r 
for that single measurement. Hence the qubit state r is conditioned on the result 
of the measurement of x. If this state is an entangled one then we have created an 
entangled state out of the unentangled initial condition l^f^f) by making the measurement. 
Furthermore, the creation of the entangled state, although probabilistic, is heralded on 
the result of the x- measurement, i.e. the result of the measurement tells us which peak 
of the distribution we are in, and hence whether or not we were successful in our attempt 
to create the state. Equations (!8!)- (fT2l) show that if the measurement result corresponds 
to the field state |$o(^)), the qubit state will be a superposition of the states \gg) and 
|ee), with a relative phase given by = 2{6 + 7r/2+u!t). Although this state is maximally 
entangled for any given phase, a two-qubit state with unknown phase has no extractable 
entanglement and can be regarded as mixed. Thus we specify a particular target state 
with a chosen phase and try to produce this given state. 

We would like to know how efficient this protocol is at producing a given target 
state. We can ask this question in the following way: if a projective measurement is 
made on x, what is the probability that the state of the qubits after the measurement 
will have a fidelity (given hj F = |(a|6)p for two pure states \a) and \b)) greater than 
Fmin with our target state? We call this value the probability of success, Pg. 



where 6(F — Fmin) is zero for F < Fmin and unity for F > Fmin- This criterion both 
allows us to specify how good our state is and tells us how often the procedure works. 
In figures m and [5] Ps is plotted for the maximally entangled state (|ee) + \gg))/^/2 as a 
function of time, with a minimum fidelity Fmin = 0.9. 

We find that Pg shows a series of fiat-topped peaks as a function of time, figure HI 
occurring at twice the resonator frequency uj, indicating that the there is a significant 
chance that the measurement will produce the state (|ee) -|- \gg))/^/2 twice every 



r 




(14) 
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Figure 4. Revival of the qubits' initial state and the probability the entanglement 
procedure is successful, P,- The black curve shows the probability that the qubits 
are in the state \gg)^ tracing over the resonator. The (lighter) red curve shows the 
probability of obtaining the state {\gg) + |ee))/\/2 with a fidelity F > Fmim Ps, as a 
function of time. The vertical dashed lines indicate the revival and sub-revival times 
at tr = 2tt^/E/X and tr/2 respectively. We have chosen X = uj = 1, n = 200 and a 
minimum fidelity Fmm = 0.9. 



oscillation cycle. Note that the probability of producing this state is quite large (close 
to 0.5) for a significant period of time. We obtain the state (|ee) + \gg))/\/2 when the 
measurement of x corresponds to the peak with frequency u (the central peak in figure 
[3]d). The peaks are most widely separated at the time tr/4, but it is worth emphasising 
that we observe peaks in P, close to 0.5 as soon as the peaks do not overlap, i.e. at a 
much shorter time than the revival. 

At other points in the oscillation cycle, the qubits will have a different relative 
phase. In figure [5] we have also plotted the probability of obtaining the states 
{\ee)+e-''^\gg))/V2with(f) = 0,7r, ±7r/2. We can therefore produce any of the maximally 
entangled states of the form {\gg) + e~^'^\ee))/ \/2 by making the projective measurement 
at the appropriate time. 

We might expect that as the phase is evolving continuously, we would only get the 
desired target state at instantaneous points in time. The fact that we can still get a non- 
zero Ps for a finite time is due to two factors: firstly that the states {\gg) + e~^'^\ee)) / 
for different 0's are not orthogonal, and secondly that the coherent state has a finite 
size. 

States with different values of are not in general orthogonal. Of course, if we 
were to require that we obtain our target state with unit fidelity, then we would expect 
to obtain our target state only at a single moment in time, and indeed the peaks in Ps 
get narrower as Fmin is increased. However, if we choose an Fmin < 1 which we consider 
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Figure 5. The probability, Ps, of achieving a fidehty greater than 0.9 for target 
states of the form {\gg) + e*"^|ee))/-\/2 with phases 0,n,±n/2. The coupHng is given 
hy X = u! ~ 1, and h — 200. The curves all reach a constant value close to 0.5 for a 
finite period. The value of exp(i^) = ±1, ±z for the state produced is indicated above 
each peak. 



"good enough," then a range of values of (f) will have fidelity greater than this with the 
target state. This leads to a simple trigonometric relationship between the width of the 
peak in Pg and the desired fidelity, 

^^oo _ arccos(2F^i^ - 1) ^^^^ 

which is shown in figure [6l In the limit n oo, we find that P, as a function of time is a 




Minimum Fideiity 



Figure 6. Full width at half maximum height for the Pg peak, as a function of fidelity. 
The numerically calculated value for n = 300 (light blue dashed line) is compared to 
the analytical form in equation (fT5|) (black line) . We get good agreement between these 
even for values of n where the peak shows significant broadening because the shape 
of the peak (figure [7]) changes in such a manner that the FWHM value is relatively 
unaffected. 



series of top-hat functions with a width given by equation (fT5l) and a height of 0.5. The 
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projective measurement "perfectly" produces a state of the form {\gg) + e~'"^|ee))/v^ 
and the width of the peaks in Pg is solely due to the finite overlap between this state 
and the target state. 

For finite n coherent states we see in figure [7| that the value of Pg changes 
continuously, rather than as a top-hat function. This smoothing is due to the fact 
that the coherent state has a finite width. If the centre of the peak in the x-distribution 
corresponds to the state {\gg) + e~*'^|ee))/v^ the leading and trailing edges of the peak 
correspond to qubit states with phases slightly above or below 0. This means that even 
when the centre of the peak has a fidelity lower than Fmm, the leading or trailing edge 
may have the target value of 0, leading to peaks in with a greater width, as can be 
seen in figure [3 The difference in phase across the distribution also means that when 
the centre of the peak has the desired phase, the edges of the peak have the "wrong" 
phase and Pg will be reduced (observable as a reduction of the step height in figure [7]) . 
As n is increased, the difference in phase between the edges and centre of the peaks 
disappears, so that for n — oo, the whole of the peak corresponds to the exact state. 
Numerical calculations indicate that the width of the peak is approximately given by 
the "ideal" width At°°, plus a term due to the finite size of the coherent state, 

^tpWQM = 7^ h At (16) 



where K{Fmin) is a constant that depends on F„ 




Figure 7. The probability, P, of achieving desired fidehty Fmin as a function of 
time over one peak, plotted for coherent states with n = 300, 200, 100, 50, 25 (top to 
bottom). The two sub-plots show Pg for different values of Fmin- For large coherent 
states and low Fmin the peaks approach a top-hat form. For smaller coherent states 
or larger Fmin the states become more rounded. 



There are some points in the evolution that do not lead to a high value of Pg for 
any phase 0. This can be observed in figure [5l where the peaks in Pg for the states 
with phases = 0, ±7r/2 overlap (indicating a continuous change from one state to the 
other), but Pg is zero for a finite period between the peaks representing states with 
(j) = IT and = ±7r/2. These periods, which occur four times a cycle, represent the 
times when two of the x-distribution peaks overlap and so the qubit state produced by 
the x-measurement is not of the form {\gg) + e~^'^\ee)) /\/2. 
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In the limit of n — > oo, the peaks in Pg are step functions of height 0.5 with a width 
determined by equation (fT5|) . which means we can achieve any desired fidehty with a 
probabihty of close to 0.5 as long as we can measure x at an exact time. Of course, we 
do not have instantaneous measurements, and so, depending on the desired fidelity and 
the time resolution of the measurement, it may be desirable to have a smaller coherent 
state. As -F^m increases, we see that the width of the step function decreases (figure 
|6]), meaning that the measurement must occur within an increasingly specific window 
in time. If this becomes a limiting factor, it may be preferable to use a smaller coherent 
state, increasing the window in which a high fidelity state can be found at the expense 
of reducing the probability of obtaining that state. An additional problem in any real 
measurement will be imprecise in space as well as time, i.e. the value of x returned by 
the measurement is probabilistic. It is clear that if the imprecision in the measurement 
is much smaller than the widths of the peaks in the x-distribution, then the imprecision 
will have little effect, and we can treat the measurement as projective. However, we 
also find that as long as the imprecision is smaller than the separation between peaks, 
the protocol still produces high fidelity states, although the effect of finite n in reducing 
the fidelity becomes more pronounced. 

It should also be noted that this technique can easily be scaled up to entangle larger 
numbers of qubits. For larger numbers of qubits the entanglement of the states created 
is not as easy to quantify, and it may not be possible to couple them to the resonator 
with equal strength, and so further research is required in this area. 

4. Conclusions 

We have described a system of two qubits capacitively coupled to a superconducting 
microwave resonator, and shown that we can devise a protocol for heralded probabilistic 
production of entangled states of the qubits. This method only requires an initial 
preparation of the system in a product state of the qubit ground states, a coherent state 
of the field in the resonator, and the ability to perform a homodyne measurement on 
the resonator. This protocol can produce states with a fidelity as high as desired with 
a probability approaching 0.5 in the limit of an infinitely large coherent state. We have 
shown that there is a trade-off between the fidelity desired, the probability of obtaining 
the desired state and the time window in which the measurement must be performed. 
Choosing a smaller coherent state increases the time window of measurement but reduces 
the probability of obtaining the desired fidelity. 

We have shown that this protocol is in some sense the dual of the phenomenon 
of revival in that both cases rely on the fact that we can regard the whole system as 
consisting of several field coherent states for each qubit state, with the different field 
states oscillating with different frequencies. Revival occurs when the time evolution of 
these states at these frequencies brings them all into phase with each other. In contrast 
the entanglement protocol works best when the states are well separated in phase. 

Being based on measurement of a quantum mode, which is coupled separately to 
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both qubits, our approach to generating entanglement contrasts with approaches where 
the two qubits, ahhough not interacting directly with each other, are both coupled 
directly to the same measurement apparatus or detector. A solid state qubit example 
of this is given by Ruskov and Korotkov |33j, for two quantum dots coupled to a 
point contact, or two Cooper-pair boxes coupled to a SET. Our entangling protocol 
also contrasts with approaches such as that of Schneider and Milburn|34]. where the 
common resonator mode is both damped and driven for all times. In such approaches 
feedback and control[35l [32], based on homodyne measurement, can be utilised to 
enhance the results. It is possible that similar application of feedback could help for our 
approach. However, even without any enhancement over our present results, it should 
be noted that such high fidelity, but probabilistic entanglement — heralded through the 
measurement outcome — has good use in Quantum Information Processing. Given the 
recent progress with superconducting qubit experiments, initial investigations of such 
entangling approaches should soon be possible. In the longer term, probabilistic but 
heralded entanglement generation can be used [Ml ET] as a basis for efficient quantum 
computation using the cluster state approach j38]. 
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